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Abstract. We have made 12 CO 2-1 and 1-0 maps of eleven molecular outflows associated with intermediate to 
high-mass young stellar objects (YSOs) in order to establish whether the correlations between outflow parameters 
and source bolometric luminosity hold in the high-mass regime. It is important to consider the effects of Malmquist- 
type biases when looking at high-mass YSOs, as they are generally much more distant than their low mass 
counterparts. We therefore chose only objects located at ~2kpc. We find that the relations show much more 
scatter than is seen in similar studies of low-mass YSOs. We also find that the mass-spectrum is significantly 
steeper in high-mass outflows, indicating a larger mass-fraction at lower velocities, a low collimation factor (~ 
1-2) and no Hubble-like relationship. 
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1. Introduction 

It has been well established that outflows represent an im- 
portant evolutionary stage in the formation of low-mass 
stars (Bachiller, 1996), but systematic studies to iden- 
tify bipolar outflows from massive YSOs have only begun 
much more recently. There is now growing evidence that 
bipolar molecular outflows are ubiquitous with YSOs of 
all masses. Shepherd & Churchwell (1996a) searched for 
high velocity (HV) CO wings toward ~ 120 high- mass 
star forming regions, and found that HV gas was present 
in ~ 90% of sources, indicating that molecular outflows 
are also a common property of high-mass YSOs. Most 
luminous YSOs seem to have outflows with much wider 
opening angles than their low-mass counterparts (Richer 
et al. 2000), but this may be due to the poor resolution 
that the known flows have been observed with, and/or the 
poor statistics available due to the small number of flows 
which have been studied. It is difficult to make predictions 
about the general properties of high-mass flows due to the 
selection effects inherent in the small existing samples. 

Shepherd & Churchwell (1996b) determined the ener- 
getics of 5 massive YSO outflows and showed that they 
fit well onto the correlation between outflow force (Fcoj 
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sometimes referred to as momentum flux) and bolomet- 
ric luminosity presented by e.g. Cabrit & Bertout (1992; 
hereafter CB). Henning et al. (2000) studied three more 
massive outflows with similar results. However, we have 
shown (Moore & Ridge, in preparation) that these corre- 
lations, which use heterogeneous, unconstrained samples, 
may be contaminated by biases due to source distance. 
We have therefore undertaken a program of observations 
to study a larger sample of high-mass outflows at a single 
distance of 2kpc, with consistent sensitivity, spatial reso- 
lution and analysis. The sample is discussed in detail in 
section 2, then the observations and results are presented 
in sections 3 and 4. Sections 5 and 6 contain the discussion 
and conclusions respectively. 

2. The Sample 

We selected a base sample of 25 objects from the com- 
pilation of Wu et al. (1996), with the selection criteria 
Lboi> 10 2 L Q , and distance 2.0±0.3kpc. This full sam- 
ple was reduced to a final list of 11 targets by rejecting 
objects known to have very weak and diffuse flows (e.g. 
AFGL2591), complex and possibly multiple flows (e.g. 
DR21), and flows with little obvious bipolarity (i.e. flows 
which may be pole-on, for which it is difficult to define 
the flow extent). Table 2 lists the selected outflows. All 
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Source 


R.A.(1950) 


Dec.(1950) 


L 
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hh:mm:ss 


dd:mm:ss 
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kms -1 
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10 b 


1.7 
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NGC 63341 


17:17:34 


-35:43:47 


8.00 


X 


10 4 


1.74 


-7.5 
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GGD27 


18:16:13 


-20:48:43 


2.0 


X 


10 4 


1.7 


+12 
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S88B 


19:44:40 


+25:05:30 


1.80 


X 


10 5 


2.0 


+21 


1 


IRAS 19550 


19:54:59 


+32:48:29 


1.50 


X 


10 2 


2.0 


+ 12 


1 


IRAS 20188 


20:18:51 


+39:28:18 


1.26 


X 


10 4 


2.0 


+2 


1 


W75N 


20:36:51 


+42:27:20 


1.40 


X 


10 5 


2.0 


+9 


1,2 


W3IRS5 


02:21:53 


+61:52:21 


1.10 


X 


10 6 


2.3° 


-40 


1 


AFGL437 


03:03:32 


+58:19:37 


2.40 


X 


10 4 


2.0 


-39 


1,3 


AFGL 5142 


05:27:28 


+33:45:37 


3.80 


X 


10 3 


1.8 


-4 


1,4 


AFGL5157 


05:34:33 


+31:57:40 


5.50 


X 


10 3 


1.8 


-18 


1 



°Imai et al. (2000) find a distance to W3IRS5 of 1.8±0.14kpc based on maser observations. 6 Refs: 1. Wu et al. (1996), 2. Moore 
(1989), 3. Weintraub & Kastner(1996), 4. Hunter et al. (1995). 

Table 1. Single distance sample of outflow sources. Columns 2 & 3 give the positions of the sources. All offsets quoted 
throughout this paper are given relative to these positions. Vlsr is the velocity of the ambient cloud material relative 
to the local standard of rest. Lbol is the bolometric luminosity of the source believed to be the driving source for 
the outflow. Luminosities given by Wu et al. (1996) were checked by referring to the original paper, and found to be 
unreliable in several cases. 



of these flows have been previously observed with low 
resolution and with different sensitivities, but few have 
been previously mapped in any detail. The sources are 
selected from an existing, heterogeneous catalogue of out- 
flow sources, and therefore may not be truly representative 
of the population properties. However, by selecting objects 
with a range of luminosities, all at the same distance, bi- 
ases should be reduced to minimum. The bolometric lu- 
minosities listed in Table 2 are generally based on IRAS 
fluxes and thought to be reliable. 

3. Observations and Data Analysis 

Observations were carried out between 1998 August and 
1999 June at the 15m James Clerk Maxwell Telescope 
(JCMT) on Mauna Kea, Hawaii and at the NRAO 12m 
Radio Telescope (NRAO) at Kitt Peak, Arizona. 

3.1. JCMT Observations 

We obtained 12 CO J=2^1 "flexibly-scheduled" observa- 
tions between 1998 August and 1999 May and 4 sched- 
uled nights of observations in 1999 July. At the observ- 
ing frequency of 230 GHz the half-power beam width of 
the telescope is 21". We used the DAS spectrometer with 
1657 channels, providing an effective velocity resolution of 
0.4kms _1 . The system temperatures during the observa- 
tions varied between 300 K and 500 K depending on the 
weather. Pointing and focus checks were carried out reg- 
ularly on planets or bright galactic sources. Pointing was 
found to vary by less than 5" during the runs. 

The observations were carried out in raster-mapping 
mode with a grid spacing of 10" producing fully sampled 



maps. In this mode, the telescope observes one "off" po- 
sition per row of on-source observations. Maps were re- 
peated 4 times to obtain the required rms sensitivity of 
0.2 K per lkms -1 resolution element. The reference po- 
sitions were checked for the presence of CO emission by 
single position-switched observations. 

The data were calibrated using the standard chopper- 
wheel technique, so the intensity scale for all JCMT 
spectral data presented here is expressed in units of T* 
(Kutner and Ulich, 1981) which is the source antenna tem- 
perature corrected for atmospheric and rearward spillover 
losses. This is converted to radiation temperature, by cor- 
recting for the forward spillover efficiency r)f ss . Single 
position-switched observations of a standard source (W3) 
were also obtained so that the data could be compared 
with the NRAO observations. 

3.2. NRAO Observations 

We obtained 13 CO J=2-l and 12 CO J=l-0 observations 
at the NRAO 12m telescope during scheduled runs in 
1998 November, 1999 March and 1999 April. These data 
were necessary in order to investigate variations in opti- 
cal depth and excitation which can significantly affect the 
masses derived via the LTE approximation. 

The telescope has a half-power beam width of 27" 
at 230 GHz and 55" at 115 GHz. At 230 GHz the dual- 
channel single-beam SIS receiver was used with two fil- 
ter bank spectrometers in parallel, each with 128 chan- 
nels. The filter banks were used at 500 kHz and 1 MHz 
bandwidths providing a spectral resolution of 0.7kms _1 
and 1.3 km s -1 respectively at 230 GHz. System tempera- 
tures varied during the runs, and channel 2 of the receiver 
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achieved consistently worse performance during the obser- 
vations. 

The observations were made in On-the-Fly mapping 
mode. As at the JCMT, the chopper wheel calibration 
method was used, yielding intensities on the T* scale. 

The mapped fields were the same size as the JCMT 
fields, plus a 30" "ramp-up" distance added to the end 
of each row. A row spacing of 8" and row scanning rate 
44"s _1 were used, resulting in fully-sampled maps. Maps 
were repeated, scanning alternately in the RA and Dec di- 
rections to remove artifacts in the scanning direction, and 
later combined and gridded together to improve signal-to- 
noise in the final map. 

3.3. Analysis 

A local thermal equilibrium (LTE) approximation was 
used to convert the 12 CO intensities into column densi- 
ties, using the 13 CO 2-1 data to correct for optical depth 
variations with velocity and position, assuming a solar- 
system isotopic abundance ratio of 89. Smaller values of 
this ratio have been observed by e.g. Langer & Penzias 
(1990) but this will affect the derived masses by less than 
a factor of two (Ridge, 2000). 12 CO 1-0 data was addition- 
ally used to calculate excitation temperatures to improve 
the analysis. 12 CO column densities were then converted 

abundance Galactic 



Xh 2 



to masses, using a standard 
ratio. 

Flow dynamics were derived following the method of 
Cabrit & Bertout (1992). A full description of the data 
reduction and analysis procedure is given in Ridge (2000). 
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AFGL 5142 
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AFGL 5157 
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IRAS 19550 
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Table 2. 12 CO flux at the positions of maximum inte- 
grated intensity of the red- and blue-shifted emission. 
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NGC 63341 
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-40.0^ 


-10.0 


NGC 6334B 


3.0^ 


6.0 


-18.0- 


-6.0 


S88B 
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AFGL 437 
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-48.0 - 
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AFGL 5142 
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AFGL 5157 


-16.0-* 


-3.0 
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W3IRS5 


-39.0^ 


-14.0 


-61.0-> 


-43.0 


GGD27 


14.5-> 


20.0 


1.5-> 


10.5 


IRAS 19550 


14.0-> 


21.0 


5.5^ 


11.5 


IRAS 20188 


4.0^ 


16.0 


-20.0- 


-2.0 


W75N 


15.0-> 


32.5 


-30.0^ 


-18.0 



Table 3. The ranges of red and blue shifted emission over 
which the integration was carried out. Velocities are LSR. 



4. Observational Results 

4.1. Morphology 

Maps of the 12 CO J=2-l integrated intensity from JCMT 
are shown in Figs. 1 and 2. Blue-shifted emission is shown 
as dashed lines, and red-shifted emission is shown as solid 
lines. Contour levels are as given in the captions, with the 
peak-flux (integrated 12 CO flux at the position of max- 
imum integrated intensity) summarised in Table 2. The 
integration limits for the red and blue shifted gas are 
summarised in Table 3. In general the flows are not well- 
collimated, but bipolar structure is evident in most cases. 
We discuss each object individually below. 

4.1.1. NGC6334I 

NGC 6334 is a large complex of H n regions and an as- 
sociated molecular cloud in the southern Galactic plane, 
containing numerous sites of recent or ongoing star forma- 
tion. NGC 63341 is the brightest far-infrared source in the 
northern part of NGC 6334 (Loughran et al., 1986) and 
is the dominant source at millimetre and sub-millimetre 
wavelengths (Sandell, 2000). The high- velocity outflow 
was first detected by Bachiller & Cernicharo (1990). 
NGC 63341 is also a rich source of molecules (Bachiller 



and Cernicharo, 1990) and has water, methanol and OH 
masers associated with it. Three ammonia masers have 
also been detected near the ends of the outflows (Kraemer 
and Jackson, 1995). 

The integrated contour map of the NGC 63341 outflow 
is shown in Fig. 1. There are two distinct flows, indicated 
by two clearly separated lobes present in the red-shifted 
emission (solid contours). It is less clear in the blue-shifted 
emission, but there seems to be an extension to the north- 
west which could be a second flow lobe. The second, less- 
dominant flow seems to be centred to the west of the 
stronger flow, and probably originates from a second ob- 
ject. This is confirmed by McCutcheon et al. (2000) who 
detect this feature in both 12 CO 3-2 and in CS, and assign 
it the name I(NW). I(N) is a third object in this region, 
located - 100" north of NGC 63341. The bottom edge of 
the outflow from this source can be seen on the north- 
ern edge of the frame. Outflow dynamics were calculated 
from the dominant outflow, closest to the IRAS source. 
Extremely high velocity emission is seen in the spectra of 
this source and there is a sharp absorption feature at Vlsr 
+6.5 km s -1 . This is a real line in the same sideband, not 
due to emission in the off position and is also seen in J=3- 
2 emission by McCutcheon et al. (2000) widespread over 
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Fig. 1. Maps of the integrated 12 CO J=2-l emission from 6 sources made at the JCMT with 21" resolution. Contours 
are at 10% intervals in all cases. Blue-shifted material is indicated by dashed contours and red-shifted material by solid 
contours. Crosses indicate the position of the infrared source associated with each of the flows. Velocity integration 
ranges are given in Table 3. 



the northern part of the NGC6334 complex, suggesting 
that a foreground cloud at a lower temperature may be 
obscuring the radiation from NGC 63341. 



4.1.2. NGC6334B 

The second source in the NGC 6334 complex, NGC 6334B 
was observed by Phillips & Mampaso (1991). Their ob- 
servations centred on the CO peak position of Dickel et 
al. (1977), and they did not find a bipolar structure in 
the high- velocity gas. However, there is an extremely lu- 
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minous FIR source (L~ 4 x 1O 5 L ; Loughran ct al. 1986) 
and H n region 40" to the north of this position, and our 
maps were therefore offset north 40" from the CO peak 
position. 

The integrated contour map of the NGC 6334B out- 
flow is shown in Fig. 1. This flow shows two lobes well 
offset to the north-east from the centre of the map. The 
peak emission of the two lobes is separated by one beam 
width. The blue lobe may stretch beyond the map bound- 
ary to the north. The flow lobes are not well collimated. 
Both the 12 CO and 13 CO spectra show peaks at and 
-5kms _1 , indicating two separate cloud components, and 
the spectrum from the red peak position shows a possible 
third peak at 3kms~ 1 . This is clearly a complex region. 

4.1.3. S88B 

S88B is a compact Hn region located within the S88 
cloud. Single-dish molecular line observations show that 
the S 88B complex is embedded within a molecular cloud 
core of relatively low mean density, n(H2) ~ 6x 10 3 cm~ 3 , 
out of which the exciting stars of the H n regions pre- 
sumably formed (e.g. Phillips & Mampaso 1991 and refs. 
therein). The outflow was mapped by Phillips & Mampaso 
(1991) with 30" resolution. 

The integrated contour map of the outflow from S88 B 
is shown in Fig. 1. Two lobes are clearly visible, with 
the red-shifted lobe (solid contours) appearing more colli- 
mated than the blue-shifted lobe (dashed contours). Two 
distinct peaks are present in the red-shifted emission, but 
the blue-shifted emission gives no indication of a second 
bipolar outflow to the south. 

4.1.4. AFGL437 

AFGL437 is a compact cluster (Kleinmann et al., 1977) 
containing at least 3 highly polarised YSOs (Dyck and 
Lonsdale, 1979). Asymmetric profiles were detected near 
this source by Schneps et al. (1978), but were not 
initially interpreted as a bipolar outflow. Arquilla & 
Goldsmith(1984) made observations at 50" resolution and 
a higher sensitivity and interpreted the profiles as pro- 
duced by a large-scale bipolar outflow, with the red and 
blue lobes located at the north and south of the central 
cluster, and their maxima separated by ~ 1.2'. Gomez 
et al. (1992) mapped the region at higher resolution and 
found a compact outflow with a very low degree of colli- 
mation. 

The integrated contour map of AFGL437 is shown 
in Fig. 1. This is a poorly collimated compact flow. 
However bipolar structure is still evident, with the blue 
lobe stretching south and the red lobe extending to the 
north. The north lobe particularly appears elongated per- 
pendicular to the outflow axis. The AFGL437 infra-red 
cluster, indicated in the diagram by the cross is not at the 
centre of the outflow. We suggest that the contour pattern 
here is consistent with the base of a very wide parabolic 



outflow in a clumpy medium. Weintraub & Kastner (1996) 
use infrared polarisation data to explain the unusual mor- 
phology of this CO outflow by the presence of an obstruc- 
tion to the north of AFGL437 which will defocus and 
deflect the flow. It may also be due to the position of the 
cluster near to the edge of the molecular cloud (Gomez 
ct al., 1992). 

4.1.5. AFGL5142 

The star-forming region AFGL 5142 contains an IRAS 
source with a bolometric luminosity 3.8xl0 3 L© 
(Carpenter et al., 1990). The CO outflow was found 
by Snell et al. (1988) with moderate resolution obser- 
vations. Its radio continuum flux density is consistent 
with a B2 ZAMS star, assuming optically thin free-free 
emission (Torrefies et al., 1992a). 

There is a cluster of infrared sources in the centre of the 
region with a radius of 0.3 pc (30" at 2 kpc; Hunter et al. 
1995), many of the sources showing a strong infrared ex- 
cess. Two cluster members, IRS 1 which has the strongest 
IR excess and IRS 2, the brightest member at K-band (2.7 
magnitudes brighter than IRS 1) have received the most 
attention in previous studies. IRS 2 is at the centre of the 
large scale outflow (Hunter ct al., 1995) and also coincides 
with the IRAS point source position (marked with a cross 
in the figure). Hunter et al. suggest that IRS 2 is in an ad- 
vanced stage of protostellar evolution, in a region where 
the molecular cloud has been dispersed. IRS 1 is thought 
to be in a much earlier evolutionary phase, coinciding with 
an unresolved thermal continuum source, a cluster of H2O 
masers and a small dense molecular cloud core. 

The integrated contour map of AFGL 5142 is shown 
in Fig. 1. This is another relatively compact flow, with 
the blue and red peaks superimposed suggesting a flow 
at a small angle to the line of sight. This flow is poorly 
collimated, and Hunter et al. (1995) suggest it may be 
multi-polar, with two flows orientated almost perpendic- 
ular to one another. Our observations do not support this 
suggestion. 

4.1.6. W3IRS5 

The source IRS 5 in the W3 giant molecular cloud was 
discovered by Wynn- Williams et al. (1972). It has a very 
high luminosity (~ 10 6 L Q ), extremely steep spectrum and 
a very deep silicate absorption feature at 9.8 fxm (Willner, 
1977). The outflow was first identified by Bally & Lada 
(1983). 

The integrated contour map of W3 IRS5 is shown 
in Fig. 1. Although relatively poorly collimated, bipolar 
structure is clearly evident in this outflow. Emission is 
present out to > 20kms~ 1 in this flow. This source shows 
strong self-absorption, especially noticeable in the spec- 
trum from the blue lobe, where the line centre, indicated 
by the peak of the 13 CO spectrum, is almost totally self- 
absorbed. 
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Fig. 2. Maps of the integrated 12 CO J=2-l emission from 5 sources made at the JCMT with 21" resolution. Contours 
are at 10% intervals in all cases. Blue-shifted material is indicated by dashed contours and red-shifted material by 
solid contours. Crosses indicate the position of the infrared source associated with each of the flows. Integration limits 
for the red and blue wings are given in Table 3. 



4.1.7. AFGL5157 

AFGL 5157, situated in the Perseus arm (Roberts, 1972), 
is an extensively studied star-forming region. Two optical 
reflection nebulae, NGC 1985 and GM 39, have been ob- 
served in AFGL 5157 (Higgs 1971, Torrelles et al. 1992a). 
Near-infrared observations by Torrelles et al. (1992b) pre- 



sented two infrared sources, IRS 1 and 2. IRS 1 appears to 
be associated with the infrared source IRAS 05345+3157, 
which has a far-infrared luminosity of 5.5 xlO 3 L© (Sncll 
et al., 1988). IRS 2 is located close to the reflection neb- 
ula GM39. In the K' imaging survey by Hodapp (1994), a 
cluster of stars can be seen in the IRS 1 position associated 
with dispersed nebulosity. 
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The IRAS source 05345+3157 is not thought to be 
the powering source of the outflow in AFGL5157, as it 
is located ~ V southwest of the outflow and molecular 
core centre (Verdes-Montenegro et al., 1989). Torrelles et 
al. (1992b) suggest that the outflow is in fact powered by 
a deeply embedded source at the centre of the ammonia 
core and associated with the strong H2O maser. 

The integrated contour map of AFGL5157 is shown 
in Fig. 2. This flow shows a high degree of collimation 
and bipolarity with weak extensions stretching for an ar- 
cminute to the east and west. The majority of the emission 
is however concentrated close to central source. There may 
be a weak second flow located to the north-west. 

4.1.8. W75N 

W75N is a well studied active star-forming region con- 
taining compact Hil regions, OH and water masers, 2 /im 
emission and about 40 embedded sources, as well as the 
well-known CO outflow (Moore et al. 1991a, Moore et al. 
1991b). It is thought to be forming several B-type stars 
(Hunter et al., 1994). The source associated with IRS 1 
and Hn region 'B', with a luminosity of ~ 1.4 x 10 5 L Q 
(Moore et al., 1991a) is thought to be driving the outflow 
(Moore et al., 1991b). H 11(B) has since been resolved into 
three subcomponents by Hunter et al. (1994). 

The integrated contour map of W75N is shown in 
Fig. 2. There is a well collimated red-shifted lobe which 
stretches south west from the source. It appears to start 
very narrow, then open up as it gets further from the 
source. 

It is more difficult to determine the morphology of the 
blue-shifted lobe, as the spectra are contaminated by emis- 
sion from the nearby cloud containing the object DR21 
which has an LSR velocity of — 3kms _1 and overlaps the 
W75 cloud over most its extent (Dickel, Dickel & Wilson 
1978). The 13 CO spectra show the two separate objects 
clearly, but it is impossible to determine how much of the 
12 CO emission is from the wings of W75 N and how much 
from DR 21. The peak at an offset of (0,-15) arcseconds is 
strongest at ~+3kms _1 and is therefore likely to be due 
to the outflow. However, this peak is compact, showing 
virtually no extension along the axis of the red outflow 
lobe. Davis et al. (1998) suggested that the small extent 
of the blue lobe is probably due to environmental asym- 
metry, with little molecular gas to the north-east of the 
source, and supported by the fact that the asymmetry of 
the flow is not reflected in H2 observations. Due to the un- 
usual nature of the blue-lobe of this flow, all calculations 
of flow properties were made by scaling from the red-lobe 
values. 

4.1.9. IRAS 19550 

IRAS 19550 was identified as a protostellar object from its 
outflow by Koo et al. (1994). This object appears to be 
forming in isolation at the tip of a large filamentary cloud. 



IRAS 19550 shows K-band emission extended along the 
east-west direction, and I and R band images show that 
the source is made up of two weak peaks shifted symmet- 
rically from the K-band peak, probably due to scattered 
stellar light (Koo et al., 1994). 

The integrated contour map of IRAS 19550 is shown 
in Fig. 2. This flow is the most compact of those we stud- 
ied. There is also some indication of a second flow located 
to the south-west of the source and at right angles to the 
known flow, particularly evident in the red (solid) con- 
tours. The source has very weak 12 CO emission. 

4.1.10. IRAS 20188 

IRAS 20188, a compact molecular cloud identified in a 
survey for molecular line emission from IRAS sources by 
Richards et al. (1987), is located in the Cygnus region and 
is not obviously associated with any known stellar cluster 
or H 11 region. The outflow was first identified by Little et 
al. (1988). 

The integrated contour map of IRAS 20188 is shown 
in Fig. 2. The morphology of this flow is difficult to de- 
cipher. There appear to be two well-collimated lobes at 
almost 90 degrees to each other, both containing at least 
two distinct peaks, which may suggest a multiple flow or 
alternatively may be a result of multiple outflow episodes. 
One explanation of this flow is that it is the base of a 
very wide parabolic flow with an asymmetric distribution 
of material. 

4.1.11. GGD27 

CCD 27 was originally catalogued as an optical candi- 
date Herbig-Haro object (Gyulbudaghian et al., 1978), but 
has now been determined to be the centre of a highly 
active star-forming region. It contains numerous near- 
infrared point sources and extensive reflection nebulos- 
ity (Yamashita et al., 1987). The deeply embedded object 
CCD 27-ILL is thought to be the driving source of the 
outflow (Aspin et al., 1994), and the nearby source IRS 2 
is a knot of dust heated by the the hot circumstellar dust 
around GGD 27-ILL. 

The integrated contour map of GGD 27 is shown in 
Fig. 2. There are two lobes well separated in velocity (blue 
lobe integrated between 1.5kms _1 and lO.Skms" 1 , red 
lobe integrated between 14.5kms _1 and 20.0kms _1 ), al- 
though they are superimposed in space, with the blue lobe 
stretching to the north, while the red lobe stretches south- 
wards. The red lobe (solid contours) appears slightly more 
collimated than the blue lobe (dashed contours). 

4.2. Outflow Properties 

4.2.1. Flow Sizes and Dynamical Timescales 

The sizes of the flow lobes were determined from the con- 
tour maps by measuring the distance in arcsec between 
the adopted centre of the flow and the 10% contour in 
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Object 


Flow Mass 


Fco 


Leo 


Flow Size 


Flow age 


Mfl ow 




M Q 


xl0~ 3 M Q kms" 1 yr" 1 


L Q 


pc 


x 10 4 yr 


M yr" 1 


AFGL5157 


672 


131 


217 


0.69 


2.49 


0.027 


AFGL5142 


1082 


10.9 


5.47 


0.80 


3.28 


0.033 


AFGL 437 


620 


235 


750 


1.36 


5.77 


0.011 


S88B 


867 


103 


222 


1.12 


4.68 


0.019 


GGD27 


407 


14.7 


15.8 


1.44 


7.63 


0.0053 


W3 IRS5 


2016 


756 


2594 


1.73 


3.60 


0.056 


NGC 6334B 


2344 


14.1 


6.2 


0.69 


2.82 


0.083 


NGC 63341 


3317 


76.6 


71.9 


0.74 


0.60 


0.55 


IRAS 20188 


2166 


2.66 


1.10 


1.67 


4.54 


0.048 


IRAS 19550 


5.77 


0.29 


0.35 


0.58 


3.67 


0.00016 


W75N 


548 


35.2 


46.4 


1.23 


1.93 


0.028 



Table 4. Flow Properties 
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Fig. 3. (a) Flow size and (b) "Dynamical timescale" versus bolometric luminosity 



the direction of the flow axis, where the flow axis was de- 
fined as the straight line joining the blue and red emission 
peaks. The values for the red and blue lobes were then 
averaged and converted to a physical size in parsec using 
the distances to the sources listed in Table 2. The dynam- 
ical timescale, td yn was then calculated by dividing the 
flow size by the intensity-weighted velocity. It should be 
noted however that the dynamical timescale may not be a 
true "age" of the flow (see e.g. Padman et al. 1997 for a 
discussion) , particularly for sources like AFGL 437 where 
there is an obstruction or AFGL 5142 where the lobes are 
overlapping. The flow sizes and dynamical timescales are 
listed in Table 4. The flows are all of the order of 1 pc 
in size. Both Rfl ow and td yn are uncorrelated (< la) with 
bolometric luminosity (Fig. 3). 

4.2.2. Flow Mass and Energy 

Table 4 lists the dynamical age, mass, CO momentum flux 
(Fco) and power (Leo) wc have derived for all the objects 
in the sample, following the method of Cabrit & Bertout 
(1992), with corrections applied for variations in optical 
depth with velocity and position in the flow (see Shepherd 



& Churchwell 1996b). A full discussion of the procedure 
is given in Ridge (2000). 

Fig. 4 shows Fco and Leo plotted against source bolo- 
metric luminosity. Fig. 4a confirms that the force required 
to drive the outflows, Fco, is approximately 100 times 
greater than the force available in radiation pressure from 
stellar photons (Lboi/c, assuming single scattering) indi- 
cated by the solid line, and therefore radiation pressure 
cannot be the sole driving mechanism for these outflows. 
A Spearman rank test gives correlation coefficients of 0.55 
between flow force and Lboi, and 0.61 between flow power 
and Lboi- Neither of these are equivalent to a significant 
correlation (2. Oct and 2.3ct respectively). There are two 
particularly noticeable outliers which fall well below the 
dashed line indicating the CB relation in Fig. 4. The 13 CO 
emission from the object NGC 6334B is complex, and may 
be optically thick even in the wings so the values of Leo 
and Fco will be lower limits. Additionally, this flow ex- 
tended outside the mapped region due to the poor in- 
formation available in the literature about the size and 
exact location of this flow, causing the mass (and there- 
fore dynamics) to be underestimated. These factors could 
easily increase the measured values of the flow dynam- 
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Fig. 4. (a) Flow force and (b) flow power vs. Lboi for the objects in the single distance sample. The solid line in (a) 
indicates the force available in stellar photons assuming single scattering. The dashed line in (b) shows the best fit 
line to the data presented by Cabrit & Bertout (1992). 



ics by up to a factor of 10, bringing it into line with the 
scatter of the other objects. There is no similar explana- 
tion for the weakness of the flow from IRAS 20188, and 
this source appears to be driving a genuinely weak out- 
flow. Additionally, there are two sources which appear to 
be driving particularly powerful outflows, AFGL 437 and 
AFGL5157. However, both these flows are within a fac- 
tor of ^5 of CB's relation, well within the uncertainties 
expected in deriving dynamical properties. 

Therefore, within the uncertainties involved in calcu- 
lating dynamical properties from outflows, only the ob- 
ject IRAS 20188 is not consistent with the previously pre- 
sented relations between outflow force and Lboi and out- 
flow power and Lboi- 

4.2.3. Mass Outflow Rate 

Fig. 5 shows the mass outflow rate, Mf — Mf/i<jyn plot- 
ted against bolometric luminosity. The mass outflow rates 
calculated for this sample are a factor of ~10 greater than 
those derived by Shepherd & Churchwell (1996b) for ob- 
jects with a similar bolometric luminosity. This discrep- 
ancy is probably due to the different techniques used for 
determining the the velocity at which the mass is domi- 
nated by cloud core emission (low-velocity cut-off in the 
integration). Masson & Chernin (1994) and Ridge (2000) 
have shown that this can significantly influence the mea- 
sured mass in a flow, and is the largest source of error 
in mass-determination. Therefore the low-velocity cut-off 
has to be selected in a consistent manner when comparing 
samples of objects (see Ridge 2000 for a full discussion of 
the method adopted here). 

There is no significant correlation between mass out- 
flow rate and bolometric luminosity - a Spearman rank 
test gives a correlation coefficient of 0.49, equivalent to a 
1.7cr correlation. This appears to be in contradiction with 
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Fig. 5. The mass outflow rate Mf versus the bolometric 
luminosity Lboi- 

the results of Shepherd & Churchwell (1996b). However 
the relation they present shows a large scatter in Mf 
among the high luminosity objects, and would probably 
not produce a significant correlation if the objects in their 
sample with Lboi< 100 L Q were disregarded. 

4.2.4. Mass- Velocity Relation 

A number of authors have observed a power-law relation 
between intensity and velocity in the high-velocity wings 
of CO (e.g. Davis et al. 1998, Lada & Fich 1996). In 
most of the objects in this sample, at least part of the 
line profile is well described by a power-law / oc V^ 1 . 
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Fig. 6. Spatially averaged intensity vs. velocity for red 
lobe of W75 N, showing a power law relation between in- 
tensity and velocity, breaking at ~ 18km s -1 . 



Object 7bl ue 7rod 



AFGL 437 


-5.13±0.20 


-2.86±0.03 


AFGL5142 


-2.40±0.07 


-2.19±0.04 


AFGL 5157 


-1.30±0.01 


-1.59±0.03 




-2.24±0.03 


-3.79±0.14 


S88B 


-9.06±0.06 


-6.66±0.09 


GGD27 


-2.60±0.04 


-3.57±0.02 






-1.55±0.12 


W75N 




-2.54±0.02 






-5.23±0.06 


W3 IRS5 


-4.99±0.02 


-2.90±0.04 




-12.52±0.02 


-4.15±0.01 


IRAS 19550 


-1.91±0.04 


-1.52±0.03 


IRAS 20188 


-1.21±0.01 


-3.74±0.06 




-5.07±0.08 




NGC 63341 


-0.98±0.03 


-0.99±0.06 






-9.23±0.4 



Table 5. The value of 7, the power-law exponent of the 
mass- velocity relation for the objects in the single-distance 
sample. Left column gives the value for the blue wing 
(Tbiue) and right column the value for the red wing (7 ro d)- 



Least squares fits to the spatially averaged spectra are 
summarised in Table 5. Where a single power law did not 
describe the emission over all velocities, two values are 
given in Table 5, the first describing the behaviour at low 
velocities (<10kms _1 ) and the second the higher-velocity 
emission (>10kms _1 ). The slope for the red wing of the 
object W75N is shown in Fig. 6 as an example. There 
is large spread in the values of both the low- and high- 
velocity 7s measured for the flows. The mass-spectrum 
slope is steeper for the high-velocity emission in all cases 
except the red wing emission of GGD 27. The blue wing 
of W3IRS5 shows particularly unusual behaviour, with 
the mass- velocity relation steepening to a slope ~ — 12 
at a velocity of ~10kms _ , then showing a small bump 
between velocities of 14 and 16kms _1 (see Fig. 7). 




log(V - V LSR ) km/s 

Fig. 7. Spatially averaged intensity vs. velocity for blue 
lobe of W3IRS5, showing a power law relation between 
intensity and velocity, at lower velocities, but a bump at 
<~16kms -1 . 



The low-velocity mass-spectrum slopes observed for 
the majority of the flows are significantly steeper than 
the often-quoted value of ~ 1.7-1.8 for low-mass YSO 
outflows (Lada and Fich, 1996). This supports the con- 
clusions of Davis et al. (1998) who also found that out- 
flows from high-mass YSOs have a steeper mass-velocity 
relation than is observed in low-mass flows. They inter- 
preted this as evidence for the more evolved nature of the 
flows they observed. We do not find a correlation between 
mass-spectrum slope and bolometric luminosity however 
(see Fig. 8a). 

The slopes derived here cannot be considered a true 
mass-velocity relation, as in most cases, the emission is 
optically thick, and therefore mass is not directly propor- 
tional to intensity. The optical depth decreases with veloc- 
ity (Ridge 2000, Shepherd & Churchwell 1996b) and there- 
fore the correction is more extreme for the lower velocity 
emission. This means that the true slope (7) of the mass- 
velocity relation for these objects would be even steeper. 
Table 5 also shows that in all but one case (GGD 27), the 
value of 7 for the red wing is smaller than 7 for the blue 
wing. It is not clear what the cause of this effect is. Optical 
depth effects would cause red-shifted self-absorption in in- 
falling material (e.g. Myers et al. 1996). 

Fig. 8 shows the slopes plotted against bolometric lu- 
minosity, dynamical timescale td yn (flow age) and wing 
extent. We do not find the clear separation between 7s 
measured for the high velocity (> lOkms -1 ) and low ve- 
locity emission that was presented by Richer et al. (2000) 
for a sample of twenty-one outflows spanning a luminos- 
ity range from ~0.3-10 6 L© compiled from the literature. 
There is some suggestion that there is an increase in the 
magnitude of 7 with luminosity, but there arc too few 
sources in the low-luminosity range to make this conclu- 
sive. The spread of 7s does seem to increase with bolo- 
metric luminosity, similar to Fig. 4 of Richer et al. (2000). 
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Fig. 8. The slope 7 of the mass spectrum M(v) plotted 
as a function of (a) source bolometric luminosity Lboi, (b) 
dynamical age td yn and (c) wing extent AV. Triangles rep- 
resent 7 for gas with low velocity (<10kms _1 ) relative to 
the source, and squares are for gas moving at more than 
lOkms -1 . Filled symbols show 7s measured for red shifted 
emission, and open symbols show 7s for blue-shifted emis- 
sion. 



They also proposed that there exists a weak relation- 
ship between outflow dynamical age and mass-spectrum 
slope in high mass YSO outflows, with the slope becoming 
steeper with age. This is not evident in Fig. 8b, although 
we only have a relatively narrow range of flow ages so 
cannot rule this possibility out completely Fig. 8c shows 
7 plotted against AV", to test whether the value of 7 is 
affected by the detectable width of the line wings. No cor- 
relation is present between either the blue- or red-shifted 
7s and AV (<C \a in all four cases). 



5. Discussion 

5.1. The Sample 

The objects in this single-distance sample were selected 
from the catalogue of Wu et al. (1996). This catalogue 
was compiled from the literature and includes YSOs iden- 
tified from their IRAS spectral energy distribution (SED), 
radio surveys for Hll regions (e.g. Wood & Churchwell 
)and radio-continuum emission, and other outflow objects 
discovered either serendipitously or as part of targeted 
observations towards known YSO candidates. The cata- 
logue cannot therefore be considered complete, and is not 
volume-limited. Wu et al. do not discuss the possible se- 
lection effects in their catalogue. Bearing this in mind, the 
sample may not be representative of the true YSO popu- 
lation. 

In particular, there may be biases in our sample asso- 
ciated with YSOs identified from their IRAS SED and/or 
existence of radio continuum. This is likely to bias the se- 
lection to older YSOs (Class I or later), as many Class 
objects are still invisible even at far-infrared wavelengths. 
Outflows found as part of targeted observations of known 
YSOs may also introduce biases, as it is often the inter- 
esting, and therefore probably unusual (strongest), objects 
that are selected for detailed multi-wavelength study. 

However, by selecting objects with a range of luminosi- 
ties with the same distance, we believe we have achieved 
a good approximation to a representative sample of the 
properties of YSO outflows. 

5.2. Correlation between Fqo and Z_boi 

We find a marginal (2.3er) correlation between outflow 
power and source luminosity (section 4.2.2). A least 
squares fit shows a relation of Fcof* Lb o i°' 62±0 ' 21 , which 
is consistent with the slope of 0.69±0.05 for the high- 
luminosity objects in the sample of CB . Our data are off- 
set to lower Leo compared to Cabrit & Bertout's data, but 
again are not inconsistent within the uncertainties. Cabrit 
& Bertout selected well-collimated flows for their study, 
and it is therefore likely to be dominated by outflows from 
Class sources (Bontemps et al., 1996). However, YSO 
classes are only defined for low-mass objects and just a 
handful of high-mass Class candidates have been identi- 
fied to date (e.g. Molinari et al. 1998). Models of high-mass 
star-formation (e.g. Stahlcr, Palla & Ho 2000 ) predict 
there may in fact be no such thing as a high-mass proto- 
star, as the stellar birthline intersects the main-sequence 
for stars of >10 M Q . It is therefore likely that the majority 
of the sources in our sample are equivalent to Class I ob- 
jects or later. Bontemps et al. (1996) find a steeper slope 
(~1) for the Class I sources in their low-intermediate mass 
(Lboi<100) sample (consistent with the slope for just the 
lower-mass (Lboi^ 100) objects in CB's sample), giving 
an almost linear relationship between Fco and Lboi, and 
indicating that all Class I sources have approximately the 
same outflow efficiency, Fco/F ra dj of ~100, where F rad is 
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Fig. 9. Outflow efficiency versus Lboi for the objects in 
the single distance sample 



highly-collimated outflows from a cluster of young stars. 
High-mass star forming regions are also generally further 
away than the well-studied low-mass star forming regions, 
and therefore the hypothesis of multiple outflows requires 
high-resolution observations in both the mid-infrared (to 
resolve the stellar cluster) and in 12 CO (to resolve the out- 
flows) in order to be tested thoroughly. The 12 CO obser- 
vations would require mm-interferometry well beyond the 
capabilities of current facilities. Alternatively, the molec- 
ular outflow may be tracing the density structure close 
to the YSO. For instance, Mellema & Frank (1997) have 
shown that aspherical bubbles can form from the inter- 
action of a central YSO wind with a toroidal circum- 
protostellar density distribution. Observations of cloud 
core morphologies (Ridge, 2000) do not seem to indicate 
a toroidal morphology however. 



5.4. Mass-Velocity Relation 



the radiative momentum flux, equal to Lboi/c. Bontemps 
et al. also find that Class sources have an outflow effi- 
ciency of ~1000, i.e. a factor of ~ 10 more efficient than 
Class Is. 

The outflows in our sample show a range of outflow 
efficiencies (Fig. 9) from ~10 to ~ few times 10 3 (i.e. 10 
to 1000 times the force due to radiation pressure assuming 
single scattering), too large for radiation pressure to be im- 
portant in driving these flows, even if multiple scattering 
is present. If outflow efficiency is an indicator of YSO class 
(and therefore age) as Bontemps et al. suggest based on 
their observational results, then the single-distance sam- 
ple must contain a mixture of Class and Class I sources. 
There is a marginal (2. la) anti-correlation between out- 
flow efficiency and bolometric luminosity in our sample. 
This may be an evolutionary effect, as massive YSOs are 
thought to evolve more rapidly, and it is therefore likely 
that on average we are seeing them at a later stage in their 
pre-main-sequence evolution. 

The outflow efficiencies agree well with other previ- 
ous works covering several decades in luminosity (e.g. 
Rodriguez et al. 1982, Bally & Lada 1983, Lada 1985, 
Levreault 1988) and with CB who found an efficiency of 
~250 for their entire sample, while for the lower luminos- 
ity sources alone the efficiency was ~1000. 

5.3. Collimation 

The flows in the single distance sample are much 
less well collimated than the molecular outflows 
from their low-mass counterparts (collimation factors, 
R-coii = R-max/Rmin, of ~l-2 in the high-mass sample as 
opposed to >6 in low-mass objects (Lada, 1985)). There 
could be several reasons which could explain the relatively 
poor collimation of these flows. One possibility is that in 
some cases of high-mass outflows we are not seeing the 
outflow from a single star, but the superposition of several 



We find an extremely steep mass-velocity relation in a 
number of the flows in the sample, indicating a high mass- 
fraction at lower velocities. This has previously been in- 
terpreted as evidence for the more evolved nature of high- 
mass outflows, as the majority of material has had time 
to slow down after its initial acceleration due to either 
prompt or turbulent entrainment (e.g. Davis et al. 1998 
). However it is also possible that high-mass outflows are 
accelerated by a different mechanism, with the majority 
of the momentum transferred to large amounts of slow- 
moving material in the flow. 

Six of the eleven objects in the sample show distinct 
breaks in the power-law between low- and high- velocity 
has been observed in some low-mass flows (e.g. 
Lada & Fich 1996). Richer et al. (2000) found that the 
slopes (7s) for the low-velocity gas are similar in sources 
of all luminosities (0< 7 <2.5), and interpreted this as 
further evidence that a common acceleration mechanism 
operates over nearly six decades in Lboi- Richer et al. also 
found a clear separation between 7s in high- and low- 
velocity gas, interpreting this as evidence for two distinct 
outflow velocity components, corresponding to a recently 
accelerated component and a slower, coasting component. 
We found a large range of 7s for both the low- and high- 
velocity gas, and therefore my data does not support these 
conclusions. 

Models of jet-driven molecular outflows (e.g. Downcs 
& Ray 1999 , Smith et al. 1997) have tried to reproduce the 
power-law of ~1.7 observed by e.g. Lada & Fich (1996) in 
low-mass flows with some success, and the hydrodynamic 
simulation by Smith et al. (1997) also predicted the change 
in slope at high velocities, due to a jet-bow shear layer con- 
sisting of molecules which survived the jet terminal shock. 
Smith et al. also predict that 7 should steepen over time, 
possibly due to the collection of a reservoir of low- velocity 
gas. Unlike Richer et al. we do not find any correlation 
between outflow age (dynamical timescale) and the mag- 
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nitude of 7, but the data contain only a narrow range of 
flow ages, so we cannot rule this possibility out completely. 

Alternatively, Downes & Ray (1999) state that in their 
model, the molecular fraction in the jet has a significant 
influence on the value of 7 their model predicts, with 7 
increasing with decreasing molecular abundance in the jet. 
If this is the case, then the high values of 7 found for 
these high-mass outflows would imply a small molecular 
fraction in the entraining jet, and the large range of 7s (~ 
1 to ~ 12) would suggest that a large range of molecular 
abundances is possible in protostellar jets. However, their 
simulation of an atomic jet, which should give an upper 
limit to the value of 7 their model can produce, gives 7 
of 3.75, much less than the observed values even for the 
low-velocity material, and it is therefore unclear whether 
the molecular abundance of the jets has any role to play 
in the slope of the mass- velocity relation. 

5.5. The Hubble Law 

The existence of a Hubble-like relation (v oc r) for molec- 
ular outflows (e.g. Lada & Fich 1996) has been chal- 
lenged by Padman et al. (1997), who state that such a 
law may apply locally at particular bow shocks, and this 
could be interpreted as a more global phenomenon in 
small or incompletely-mapped sources. This conclusion is 
supported by the numerical simulations presented by e.g. 
Downes & Ray (1999) which show that the Hubble law 
is almost certainly a local effect in the vicinity of a bow 
shock. If a Hubble law was more global, it would be dif- 
ficult to reconcile with jet-driven models which accelerate 
the molecular material in situ rather than it originating 
at the source. The observations presented here show no 
evidence for a Hubble law in any of the high-mass ob- 
jects studied (see also the channel maps presented in Ridge 
2000). 

Models such as the simulations of Smith ct al. (1997) 
reproduce Hubble-like position-velocity diagrams well. 
However, the lack of observational evidence for a global 
Hubble-like flow puts this phenomenon in doubt. 

6. Conclusions 

Using a sample of known outflow objects at a single dis- 
tance, we have investigated the physical and dynamical 
properties of intermediate and high-mass YSOs. Our main 
results are: 

— We find a similar relation between Fqq and Lboi as 
previous studies (CB , Shepherd & Churchwell 1996b) 
but with a large scatter. 

— We find a much steeper mass-velocity relation than 
has been observed in low-mass outflows. There is also 
a large range of values of 7, from ~ —1 to ~ —12. 

— None of the objects studied show a "Hubble-law" like 
relation. 

We have found that several of the often-quoted obser- 
vational results for low-mass YSO outflows are not appli- 
cable to high-mass stars. To attempt to understand fully 



high mass star formation and their outflows, more exten- 
sive statistical studies need to be undertaken, as we have 
shown that the extrapolation of outflow properties from 
low-mass stars do not necessarily agree with observations 
of high mass stars. 
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